OPTIMAL SEQUENTIAL CHANGE-DETECTION FOR 
FRACTIONAL DIFFUSION-TYPE PROCESSES 



ALEXANDRA CHRONOPOULOU AND GEORGIOS FELLOURIS 

Abstract. We consider the problem of detecting an abrupt change 
in the distribution of a sequentially observed stochastic process. We 
establish the optimality of the CUSUM test with respect to a modi- 
fied version of Lorden's criterion for arbitrary processes with continu- 
ous paths and apply this general result to the special case of fractional 
diffusion-type processes. As a by-product, we show that the CUSUM 
test optimizes Lorden's original criterion when a fractional Brownian 
motion with Hurst index H adopts a polynomial drift term with expo- 
nent H + 1/2 after the change. 



1. Introduction 

The quick detection of an abrupt change in the behavior of a stochastic 
system is an important problem in many apphcation areas, such as quahty 
control, target tracking, navigation, seismology, bio-surveillance, computer 
security. More specifically, the problem is to find a detection rule that 
raises an alarm, as soon as possible after the change has occurred, based on 
sequential observations of the system. Thus, a good detection rule should 
have small detection delay, but also a low frequency of false alarms in its 
repeated applications. 
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There are three main formulations for the change-detection problem that 
balance the trade-ofF between these two antithetic goals. A Bayesian ap- 
proach, developed by Shiryaev [29], where the change-point is modeled as 
a random variable, and two minimax approaches, due to Pollak [23] and 
Lorden [TJ], where the change-point is considered to be an unknown, de- 
terministic parameter. For a comparison of these formulations we refer to 
the paper of Moustakides [19] and for exhaustive treatments of sequential 
change-detection to the books of Basseville and Nikiforov [1] and Hadjiliadis 
and Poor [10] . 

In the present work, we focus on Lorden's approach, which has a deep 
connection with the so-called Cumulative Sums (CUSUM) test, a detection 
structure that was proposed by Page [22] and has been very popular in 
applications since then (see for example [8]). Lorden [H] quantified the 
performance of a detection rule with its worst (with respect to the time of 
the change) conditional expected detection delay given the worst possible 
scenario until the time of the change and suggested the minimization of 
this criterion subject to an upper bound on the rate of false alarms. In the 
case of independent and identically distributed observations before and after 
the change, Lorden showed that this criterion is attained asymptotically by 
the CUSUM test. Moustakides [T7] -and later Ritov [27]- proved that the 
CUSUM test is an exact solution to Lorden's optimization problem. This 
exact optimality property was extended in continuous-time in the case of a 
standard Brownian motion that adopts a linear drift (Shiryaev [30] . Beibel 
[2]), as well as in the case of a diffusion-type process with constant "signal 
to noise ratio" (Moustakides [H]). The latter was a by-product of the main 
result obtained in the same paper, according to which the CUSUM procedure 
optimizes a modified version of Lorden's criterion in the much more general 
framework of diffusion-type processes that satisfy a "full-energy" condition. 
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In the present work, we extend the optimahty of CUSUM (with respect 
to the criterion introduced in [18]) to arbitrary continuous-path processes 
and apply this general result to the case of fractional processes of diffusion- 
type. This is a class of continuous-path processes that do not admit a 
semimartingale decomposition, thus they are clearly not in the scope of 
existing optimality results. We show in particular that the CUSUM test 
is optimal (in the above sense) for detecting a change from a fractional 
Brownian motion (fBm) to a fractional Ornstein-Uhlenbeck process (fOU), 
as well as for detecting the emergence of a linear drift in a fBm for any value 
of the Hurst index. 

Finally, we establish the optimality of the CUSUM test with respect to 
Lorden's original criterion when the observed process is a fBm with Hurst 
index H that adopts a polynomial drift term with exponent H + 1/2 after 
the change. This is one of the rare instances for which the exact solution to 
this problem is known. Of course, in the special case H = 1/2, we recover 
the optimality of CUSUM for the detection of a linear drift in a standard 
Brownian motion. 

Apart from a theoretical point of view, these extensions are also interest- 
ing due the increasing presence of fBm in a wide range of applications. This 
can be explained by the fact that fBm is well-suited to model phenomena 
that are characterized by self-similarity and/or long-memory. As a result, 
it has been used as the basic building block for models in a variety of fields, 
such as hydrology, traffic networks, finance and economics (see for example 

m, a, 0, m)- 

These diverse applications have also triggered a great interest in the sta- 
tistical inference for processes related to fBm. Indeed, Kleptsyna and Le 
Breton in [12] studied the maximum likelihood estimator (MLE) for the 
parameter in the drift of the fractional Ornstein-Uhlenbeck process with 
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Hurst index H > 1/2. Tudor and Viens [H] used Malliavin calculus tech- 
niques to study the MLE for any H £ (0, 1) and a more general class of 
fractional diffusions, in which the drift coefficient is linear with respect to 
the unknown parameter. For the same class of processes, Rao [25] stud- 
ied a sequential version of the MLE. The problem of sequential testing for 
fractional diffusion-type processes was considered by Rao in [24]. 

In what follows, we establish the optimality of the CUSUM test for ar- 
bitrary continuous-path stochastic processes in Section 2, we focus on frac- 
tional diffusion-type processes in Section 3 and we conclude in Section 4. 

2. CUSUM OPTIMALITY FOR CONTINUOUS-PATH PROCESSES 

2.1. Problem formulation. Let {^,J-) be the canonical space of continu- 
ous functions on [0,oo) that vanish at 0. We denote by {S^t} the coordinate 
process on this space and by {J-t} its natural filtration, thus ^t(<jj) := uj{t) for 
every uj £ il. and J-f ■= (t(^s : < s < t) for every t > 0, whereas '■= and 
-Fq := {0,f^}- Let Pq, Pqo be two completely specified probability measures 
on (ri, We ajSsuine tlicit Pq B^nd Poo 

are locally equivalent, i.e. mutually 
absolutely continuous when they are restricted to the c-algebra J-t for any 
< t < oo, and we denote by ut the corresponding log-likelihood ratio 



We assume that the distribution of {^j}, which we denote by P,-, changes 
at some unknown, deterministic time r G [0, oo] from Pq to Pqo- Thus, 
Pr coincides with Pqo on J^t for any t G [0, r] and with Pq on J^t ^or any 
t G (r, oo). Then, it is clear that P^ is mutually absolutely continuous with 
Poo on J^t for any t > t and that 



Ut := log 



dPo 



t>0; no := 1. 



dPoo J't 



Ut-Ur = log 



dP^ 



t > T. 



dPoo 
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The goal is to find a sequential detection rule, that is an {J^(}-stopping time 
T, that minimizes the detection delay for any change point r E [0, oo). Since 
it is not possible to do so for every r, we follow a minimax approach and 
consider the following constrained optimization problem: 

inf Jm[T] when ^ Eoo[(m)t] > 7) where 

1 + 
Jm[T] := sup esssup - E^[((u)t - {u)r) \ i>], 

T>0 ^ 

where {u)t is the quadratic variation of ut- When the observed process is of 
diffusion- type before and after the change, this problem coincides with the 
criterion introduced by Moustakides in [18] (and this explains the inclusion 
of the redundant 1/2 in its definition). When {u)t is proportional to t, it is 
equivalent to Lorden's criterion, 

inf j7L[r] when EooiT] > 7, where 

(2) 

Jl[T] := sup esssup E^[(r - Ty\ JV], 

T>0 

In both formulations, the goal is to find a detection rule that minimizes 
the worst (with respect to r) conditional expected detection delay given the 
worst possible history of observations up to the time of the change, subject 
to a lower bound on the period of false alarms. Their difference is in the way 
"detection delay" and "period of false alarms" are measured. This is done 
in terms of the actual time in ([2]) and in terms of the accumulated quadratic 
variation of the log-likelihood ratio process in ([1]). 

The latter approach also has an appealing interpretation in terms of 
Kullback-Leibler information, which becomes clear with ^ and ^ and jus- 
tifies calling (as in [18]) J^m[T] the K-L detection divergence and ^Eoo[('u)r] 
the K-L false alarm divergence of the detection rule T. However, the main 
advantage of this formulation is that its solution is known for a much larger 
class of dynamics, as we discuss below. For simplicity, in what follows we 



6 ALEXANDRA CHRONOPOULOU AND GEORGIOS FELLOURIS 

will say that a detection rule is Ju -optimal if it solves the problem defined 
in ([1]) and Ji-optimal if it solves the problem defined in ([2]). 

2.2. Main result. Let us first define the CUSUM stopping time 

Sc := inf{t > : 2/t > c}, yt := Ut - inf u^, t > 0, 

0<s<t 

where the threshold c is assumed to satisfy the corresponding false alarm 
constraint with equality, that is Eoo[(u)sc] = 7 for the problem defined in 
([1]) and Eoo [Sc\ = 7 for the problem defined in ([2]) . 

As we discussed in the Introduction, it is known that this detection rule 
is i7A/~optimal for diffusion- type processes that satisfy a "full-energy" con- 
dition and additionally ^Tl "Optimal for diffusion-type processes that have 
a constant "signal to noise ratio". With the following theorem we extend 
these optimality properties. 

Theorem 2.1. The CUSUM test is J'M-optimal, if the following condition 
is satisfied 

lim {u)t = oo Pq, Poo — 0,-S- (3) 

When in particular {u)t is proportional to t, the CUSUM test is also Jl- 
optimal. 

In order to prove this theorem, we start with the following lemma, which 
reveals the structure of the log-likelihood ratio process {ut}. 

Lemma 2.2. There exist continuous processes {Xt} and {Xt}, which are 
local martingales (vanishing at 0) with respect to Poo o-nd Pq respectively and 
have the same quadratic variation, so that 

ut = Xt-^ {u)t = -Xt + ^ {u)t, t > 0. 
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Proof. Since the likelihood ratios {e"'} and {e~"*} are continuous martin- 
gales with respect to Poo and Pq respectively, it is well-known (see for 
example, Proposition (1.6), pg. 328 in [26]) that there exist two unique, 
continuous, local martingales {Xt} and {Xt} with respect to Poo and Pq 
respectively, so that 

e''* = exp{Xi - (1/2) {X)t}, e"''* = exp{Xj - (1/2) {X)t}, 

where Xq = Xq = 0, since uq = 1. Then, taking logarithms we obtain 



ut = Xt-l {X)t = -Xt + l {X)t 



and consequently {u)t = {X)t = {X)t, which completes the proof. 

□ 



This lemma has some important ramifications. First of all, from an ap- 
plication of Wald's identity it follows that for any stopping time T with 
Ei[{u)T] < oo for ? = 0,00, 

dPoo 



1 



Eoo [("«.) 



-UT 



log 



dPr 



(4) 



and that on the event {T > r} 



- E^[(n)T - (u)t I = Er [ut - Ur) l{T>r} I 

dPr 



Er 



log 



dPr 



T 

J T 



(5) 



which highlights the connection of the performance measure Ju with the 
notion of Kullback-Leibler information. However, the most important con- 
sequence of the previous lemma is that it allows us to obtain closed-form 
expressions for the performance characteristics of the CUSUM rule and es- 
tablish its optimality, following the methodology developed in [18]. 
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Lemma 2.3. Suppose that condition holds. For any stopping time T, 
c > and r G [0, oo), 

P,(Tc < oo) = 1, Poo(r, < oo) = 1, (6) 

where Tc := T A Sc = min(r, Sc)- Moreover, on the event {Tc > r} we have 

(7) 

where the functions g and h are defined as follows: 

g{x) := e'"" + X - 1, h{x) := e"" - x - 1 , x>0. (8) 
Proof. Let us introduce the following notation: 

■.= TcAmf{t>0 : {u)t>n}, n G N. 
From an application of Ito's rule we have: 



giVT^) - giVr) = J g'{ys)dys + - J g"{ys)d{y)s- 
From Lemma 12.21 and the definition of the CUSUM statistic {yt} we have 

yt = Ut-mt = -Xt + \{u)t-mt, mt := inf Ug, 

Z 0<s<t 

and consequently {y)t = {u)t = {X)t, therefore we can write 



g{yT^) - giyr) = J g'{ys) y-dXs ^ — ^ - dmgj + - y g"{ys)d{u)s. 

Now, with a re-arrangment in its right-hand side and using the fact that the 
measure drUg is carried by the set {y^ = 0}, the previous relationship takes 
the following form: 

g{yn^) - gM = 2 I (9' + 9"){ys)d{u)s - / ' g'iys)dXs - / ' g'{0)dms 
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From ([8]) it is clear that g'{x) + g"{x) = 1 for every x > and ^'(O) = 0, 
thus 



5(2/T") - givr) 



1 



rj-in 



d{u) 



g'{ys)dXs. 



Taking conditional expectation under given Jv, on the event {Tc > r} 
we have: 



d{u) 



rpn 



g'{ys)dXs 



(9) 



Since < < c on {r < s < T" < Tc} and g' is an increasing function on 
[0,oo), we have g'{ys) < g'{c). Moreover, by the definition of the stopping 
time T" it is clear that (X)rn = (m)^™ < n, therefore 

" {g'iys)fd{X)s < {g\c)f {X)Tn < {g'{c)fn< oo. 

Then, since P,- coincides with Pq on Jv and X is local martingale (that starts 
from 0) under Pq, the second term in the right hand side of ([9]) vanishes and 
consequently we can write (on the event {Tc > r}): 



Er[ff(yr,") -5'(yr)|-7v] = E^ 



d{u) 



(10) 



Since g is an increasing function on [0, oo) and yj-n E [0, c], it is clear that the 
left-hand side of this equality is bounded by g{c) — g{0) = g{c). Moreover, 
due to condition T" converges to Tc Po-a.s. as n — )• oo. Therefore, 
letting n — )• oo in the right-hand side of (llOp and applying the Monotone 
Convergence Theorem, which we can do since {u)t is increasing, we obtain 



g{c) > E, 



d{u) 



1 

2 l{T,=oo} 



d{u) 



(11) 



From condition ([3]) and Lemma 12.21 it follows that Pq{J^ d{u)s = oo) = 1 
and consequently Pt(/^°° d{u)s = oo) = 1 since by definition Pt- coincides 
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with Pq on Tt for t > t. Then, from the law of iterated expectation we have 



PrI 



d{u) 



oo 



which imphes Pr(/^°° d{u)s = oo| Jv) = 1 and consequently Pr{Tc = oo) = 0, 
otherwise the right-hand side in (jlip becomes infinite, due to condition ([3|), 
which leads to a contradiction. 

Now, going back to ()10p and letting n — )■ oo in both sides simultaneously, 
from an application of Bounded (Monotone) Convergence Theorem for the 
left (right)~hand side we obtain the first relationship in ([8]). The second 
relationship in ([8]) and the fact that Poo (5*0 = oo) = can be shown in a 
similar way, which completes the proof. 

□ 



When we set T = in the previous lemma, we obtain some very inter- 
esting information regarding the behavior of the CUSUM test. First of all, 
from ([6|) it follows that Sc terminates almost surely for any change-point, 
r, as well as when the change never occurs. Therefore, it is almost sure that 
the CUSUM test will raise a false alarm, and this explains why a constraint 
is imposed on the rate of false alarms, not the probability that a false alarm 
will occur. 

Furthermore, from the first relationship in ([7]) it follows that on the event 
{Sc > t} we have 

since ys^ = c, due to the path-continuity of the CUSUM process, {yt}- 
Therefore, since G [0, c] on {Sc > r} and g is increasing function on 



OPTIMAL SEQUENTIAL CHANGE-DETECTION 11 

[0,00), it follows that 



Jm[Sc] = sup eSSSUp ]- £r[{{u)Sc - {u)r) l{5c>r} I ^r] 
r>0 ^ 

= sup esssup ^r[[g{c) - g{yr)) 1{&>t} I = g{c) 

r>0 



(12) 



In other words, the worst-case scenario for the CUSUM test occurs when 
the change occurs at time r = 0, i.e. J^m[T] = \ Eo[(^)t]- 

Finally, setting r = in the second relationship in ^ we obtain 

^ Eoo[(n)5j = Eoo[/i(yTj] = h{c). 

Therefore, for the CUSUM test to satisfy the false alarm constraint with 
equality, its threshold c should be chosen as the (unique) solution to the 
non-linear equation h{c) = 7. The following lemma shows that, without 
loss of generality, we can actually restrict ourselves to stopping times that 
satisfy the false alarm constraint with equality. 

Lemma 2.4. Suppose that condition ^ holds. For any stopping time T 
such that i Eoo[{u)t\ > J, there exists some c > so that ^ Eoo[{u)tJ = 7 
andJM[Tc] < Jm[T]. 

Proof. Consider an arbitrary stopping time with ^ Eoo[{u)t] > 7 and define 
the function ip{c) := ^Eoo[{u)tc]- From the second relationship in ([7]) it 
follows that il^{c) = Eoo[h{yTc)]: which implies that il^{c) is a continuous 
function. Then, since ip{0) = 0, tlj{oo) = | Eco[{u)t] > 7, there exists some 
c > so that ip{c) = 7. By definition, < T, which implies J'm[T(.] < 
J'm[T] and completes the proof. 

□ 

Proof of Theorem \2.1\ The proof is based on the fact that for any stopping 
time T and any c > we have the following lower bound 
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which can be shown in exactly the same way as in Theorem 2 of [18]. From 
(jl2p and the fact that ys^ = c almost surely it follows that both sides of the 
inequality become equal to g{c) when we set T = Sc- Due to this observation 
and Lemma 12.41 it suffices to show that 



5(c)=mf— ^-p^ when -E^[{u)r]=r 

This can be done in exactly the same way as in Theorem 3 of [IS]. □ 

2.3. The case of diffusion-type processes. Let us now illustrate Theo- 
rem [2T] in the case that the observed process {^t} is a standard Brownian 
motion before the change and adopts a (random, in general) drift after the 
change. More specifically, let Poo be the Wiener measure and suppose that 
the post-change measure Pq is induced by the following dynamics 

^t = Wt+ [ f^sds, t>0 
Jo 

where {Wt} is a standard Brownian motion under Pq and {nt} an {^t}~ 
adapted process that satisfies 

ft 



oo 



fiids <ooj =1, t> 0, (13) 
'0 ^ JQ 



exp{^ ^xs 'd^s-\ l^l ds}] =1, t > 0. (14) 



These two conditions guarantee that Pq is indeed well-defined and locally 
equivalent with Pqo- Moreover, from Girsanov's theorem we obtain the 
following explicit representation for their log-likelihood ratio, 

ut= I Hsd^s-l; I f4 ds, t > 0. 
Jo ^ Jo 

It is then clear that {u)t = /q A ds and Theorem 12.11 implies that the 
CUSUM test is j7M-optimal if the following condition is also satisfied 

oo 

fi^ds = oo, Pq, Poo — a.s. (15) 





OPTIMAL SEQUENTIAL CHANGE-DETECTION 13 

This is exactly the result obtained in |18j . When the process {fj^t} reduces 
to a constant, {u)t is proportional to t and Theorem 12.11 implies that the 
CUSUM test is -optimal, which is also a well-known fact. 

However, Theorem 12. II is a non-trivial generalization of existing optimal- 
ity results, as there are many continuous-path processes that are not of 
diffusion-type, not even semimartingales. We illustrate this point in the 
next section. 

3. The case of fractional diffusion-type processes 

Before we consider the application of Theorem 12.11 in the context of frac- 
tional diffusion-type processes, let us first define fractional Brownian mo- 
tion, the basic building block of these processes, and present its main prop- 
erties. 

3.1. Fractional Brownian Motion: A Quick Review. Let {B^^} be 
a stochastic process defined on the canonical space {^l,J-) and let P be a 
probability measure on this space. We will say that {B^^} is a fractional 
Brownian motion (fBm) under P, or simply a P-fBm, if it is a centered, 
continuous, Gaussian process with covariance structure 

^[B^B^] = l{t'^ + s'^-\t-sn, t,s>0, 

where E refers to expectation with respect to P. As a consequence of its 
definition, the process {B^^} has a number of interesting properties. Thus, 

(i) it is i?-self-similar, in the sense that {Bl^}t>o has the same finite- 
dimensional distributions as {c~^ ^rf loo ^'^^ every c > 0, 

(ii) it has stationary increments, which however are independent only 
when H = 1/2. When H < (>)l/2, they are negatively (positively) 
correlated and the process exhibits short (long) -range dependence, in 
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the sense that 

oo 

J]E[(Sf -5f„i)i?f]<(=)oo, 

n.=l 

(iii) it has Holder continuous paths of order H — e for every < e < H , 

(iv) it has finite variation (in an sense), equal to c'^t on any finite 
interval [0,t], where 

c'h:= [ \x\^/^^{x)dx, 0(x) := ^e-^'/2_ 
Jr V 27r 

(v) is not a semimartingale, i.e. it does not admit a decomposition as 
the sum of a local martingale and a term of finite variation (see, for 
example, [2S])- However, the transformed process 

Mf:=/ kH{t,s)dB^, t>0 (16) 
Jo 

is a square-integrable martingale with quadratic variation {M^)t = 
i2-2H^ where 

kH{t,s) := CjI s^~^ {t- s)^~^ , 0<s<t, 
ch and Xh are positive constants defined as follows 

CH := 2H r (3/2 - H) T (F + 1/2) , 

_2HT{3-2H)r{H + l/2) 
^ ''~ T{3/2-H) ' 

and r{x) := s^~^e~* ds is the Gamma function. This result 

was shown by Molchan [16] and more recently by Norros et. al 

pO] . where the process was called the fundamental martingale 

associated with fBm. 

For an exhaustive treatment on the properties of fBm we refer to Chapter 
5 of Nualart |21j . Here, we will only add an extension of Levy's classical 
characterization theorem to fBm that was recently established by Hu et al. 
[9] (see also Mishura and Valkeila [E]). This result essentially says that 
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properties (iii), (iv) and (v) characterize fBm. More specifically, if y is a 
continuous, centered, square- integr able stochastic process with 

(a) Holder continuous paths of order H — e for any e > 0, 

(b) finite 1 /i?-quadratic variation (in an L^-sense) that is equal to c'jjt 
on any interval [0, t], 

(c) and the process M, defined as in (|16|) with replaced by Y, 
is a martingale, whose quadratic variation is absolutely continuous 
almost surely with respect to the Lebesgue measure when H > 1/2, 

then y is a fBm with Hurst index H. 

3.2. The case of a fBm adopting a random drift. We now return to 
the setup of our change detection problem and we focus on the special case 
that the observed process {^t} is a fBm that adopts a random drift after the 
change. More specifically, in what follows we assume that {S^t} is a fBm with 
Hurst index H under Poo, whereas Pq is induced by the following dynamics 



where {fJ-t} is an {J"f}-adapted process and {B^^} a fBm under Pq with the 
same Hurst index, H. 

When H = 1/2, we recover the context of Subsection \'2.3\ where we saw 
that Pq is indeed well-defined and locally equivalent to Poo when conditions 
([T3|) and (dH) are satisfied and that the CUSUM test is J'a/ -optimal if addi- 
tionally condition (llSp holds. Our goal in this section is to obtain analogous 
conditions when H ^ 1/2, in which case {S^t} is not a semimartingale and 
the classical Girsanov theorem does not apply, thus the previous conditions 
are no longer appropriate. 

In order to do so, we will work with the transformed process 




(17) 




(18) 
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Since {^(} is a fBM under Pqo, from property (v) it follows that {Q} is 
a square-integrable martingale with quadratic variation {()t = Xjj^ t^~'^^ 
under Pqo- This property allows us to obtain an explicit representation for 
the log-likelihood ratio process {ut} using a version of Girsanov's theorem, 
although the observed process {^j} itself is not a semimartingale. This is a 
well-known result in the Stochastic Calculus with respect to fBm (see [7], 
|13j). Here, we will provide a proof that relies on the characterization of 
fBm that we previously discussed. 

In what follows, we assume that the paths of {/it} are sufficiently smooth 
so that the process 



d /■* 

Qt-=TF7r kHit,s) Hsds, t>0 
"vC/t Jo 



(19) 



is well-defined, where the derivative is understood in the sense of absolute 
continuity. 

Theorem 3.1. // the following conditions hold, 



Poo( / Qfd(C). <ooj = 1, t>0, 
ft 1 ft 







exp{^ Q,dCs-lJ^ QUiOs} 



1, t>0, 



then Pq is locally equivalent to Poo CLnd their log-likelihood ratio process ad- 
mits the following representation 

Ut = f Qs dCs - \ f Ql d{C)s, t > 0. (20) 
Jo ^ Jo 

Then, the CUSUM test is "Optimal, if additionally the following condi- 
tion is satisfied 

POO 

/ QU{C)s = oo Po, Poo -a.s. (21) 
Jo 
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Proof. Due to Theorem 12. !( we only need to show that the log-likehhood 
ratio {uf} admits representation ([20]) . In order to do so, it suffices to show 
that if the post-change measure Pq is defined by 

dPo 



dPoo 

then the process 



exp{^ Qs Ql d{Os], t > 0, (22) 

■.= it- [ l^sds, t>0, 
Jo 



is a fBm under Pq. Moreover, due to the characterization theorem that we 
discussed in the end of the previous subsection, it suffices to show that {B^} 
satisfies properties (a), (b) and (c). Since {B^^} is a "shifted" version of 
{£,t}, which is a fBm under Pqo, it clearly satisfies (a) and (b). It remains to 
show that the process { Jq knit, s) dB^} is a Po-martingale whose quadratic 
variation is absolutely continuous with respect to the Lebesgue measure 
when H > 1/2. Indeed, 

/ kH{t,s)dB^= I kH{t,s){d^s- f^sds) 
Jo Jo ^ ^23) 

= Ct- f QsdiOs 



Jo 

where the second equality follows from the definitions of ^ and Q in (|18p 
and ()19p , respectively. But from (j22p and Girsanov's theorem we know that 
the process in the right-hand side of (j23p is a Pg-martingale with quadratic 
variation {Qt = t^~'^^ , which completes the proof. □ 

In what follows, we apply this result in some interesting special cases. 

3.2.1. Fractional diffusions. The conditions of Theorem 13.11 are satisfied for 
a large class of fractional diffusions and any Hurst index H. This is the 
content of the following corollary, which is based on the results of Tudor 
and Viens [31], Lemma 3. 
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Corollary 3.2. The CUSUM is JM^optimal when fit = b{^t) and b is a real 
function so that 

(i) xb{x) has a constant sign for all x > and a constant sign for all 
x<0, 

(ii) \b{x)/x\ = c + r(x) for all x, where r{x) — )■ oo as x ^ oo. 

These two conditions are satisfied when for example b{x) = cq + cix + 
{\x\ A 1)" with Co, ci E M and a G [0, 1). When in particular a = 0, b{x) is an 
affine function and, under the post-change measure Pq, {^t} is a fractional 
Ornstein-Uhlenbeck process. This is the fractional analogue of the classical 
Ornstein-Uhlenbeck process and has been used in financial modeling (see 
for example [4] , [5] ) . For a study of its main properties we refer to Cheridito 
et al. [6]. 

3.2.2. Polynomial drift. When /Uf = t", where a is a real constant, the post- 
change dynamics (fTTl) reduce to 



ft +a+i 

6 = + / = 5f + —— 

Jo a + 1 



'0 

and the process Qt takes the following form 

Qt = -:T^ f s'^^^^-it-sf^-^ds. 

CH d{Qt Jo 

After some algebraic manipulations we obtain Qt = dH,at°', where 

_ r(3 - 2H) r(3/2 -H + a) 2 - 2H + a 
r{3-2H + a)r{3/2- H) 2 - 2H ' 

and consequently 

f Ql diOs = VH,a VH,a := , ' " ^ • (24) 

Jo Xh I- H + a 

We can now state two interesting corollaries. 

Corollary 3.3. // = t" with a + 1 > H, then the CUSUM test is Jm- 
optimal. 
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Proof. From (I24p it is clear that condition (12 ip is satisfied when a > H — 
thus the corollary follows directly from Theorem 13. 1[ 
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1, 

□ 



Corollary 3.4. If fit = with a + 1 = H + 1/2, then the CUSUM test is 
Jl -optimal. 

Proof. From (I24p it is clear that when a = H — 1/2, {u)t is proportional to 
t, thus the result follows from the previous corollary and Theorem 12.11 □ 

One of the implications of Corollarv l3.3l is that the CUSUM test is always 
c7M~optimal when a linear drift (a = 0) emerges in a fractional Brownian 
motion, no matter what the value of the Hurst index, H, is. 

Corollary 13.41 provides a class of processes for which the CUSUM test 
optimizes Lorden's original criterion. To our knowledge, such an optimality 
property has been established only in the case of diffusion-type process for 
which the quadratic variation {u)t is proportional to t (see the discussion in 
pg. 313 of [IB]). Of course, the case of a linear drift emerging in a standard 
Brownian motion corresponds to the special case H = 1/2 of Corollarv 13.41 

3.3. Extensions. It is possible to generalize Theorem 13.11 to the case that 
the pre-change measure Poo is induced by the following dynamics 



where B is a Poo^fBm and a : R-|- — t- M+ is a deterministic, non-vanishing 
real function with 5-Holder continuous paths for some 5 > 1 — H. This 
smoothness condition guarantees that the integral in (j25p can be defined in 
a Riemann-Stieljes (or Young) sense (see [32]). Then, the proof of Theorem 
13.11 goes through as long as we modify the definitions of C, and Q in ()18p and 
([19]) respectively as follows 




(25) 
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Note that when cr is a stochastic process, the integral in ([25]) cannot be 
defined as a stochastic integral in an Ito sense, since fBm is not a semi- 
martingale. 

4. Conclusions 

In this work, we extended the optimality properties of the CUSUM test 
in two directions. 

First, with respect to a modified version of Lorden's criterion, where de- 
tection delay is measured in terms of the quadratic variation (that is accu- 
mulated between the time of the change and the time of stopping) of the 
log-likelihood ratio between the post and pre-change distribution. With 
respect to this criterion, the optimality of CUSUM was established for ar- 
bitrary continuous-path processes, generalizing existing optimality results 
that refer to processes of diffusion-type. As an application, we obtained 
sufficient conditions for the optimality of the CUSUM procedure when the 
observed process is a fractional Brownian motion that adopts a random 
drift. We saw that these conditions are satisfied in some interesting special 
cases, such as when a fractional Brownian motion turns into a fractional 
Ornstein-Uhlenbeck process or it adopts a linear drift, for any value of the 
Hurst index. 

Second, with respect to Lorden's original criterion, we proved that the 
CUSUM test is optimal when a fractional Brownian motion with Hurst index 
H adopts a deterministic, polynomial drift term with exponent H +1/2. In 
this way, we generalized the well-known optimality of CUSUM in the case 
that a linear drift emerges in a standard {H = 1/2) Brownian motion. 
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